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1 Introduction 

In this paper we consider quantum cosmological models for Bianchi type I models in 
the generalized vacuum scalar-tensor theories of gravity. Interest in these theories 
has been widespread in recent years in connection with inflation and string theories. 
They are defined by the action 



A 



j d^x./^e-^[R-uj{<^){V^f ~2\{^)\, (1) 



where R is the Ricci curvature of the space-time and g is the determinant of the metric 
g^y [|l|. The dilaton field $ plays the role of a time-varying gravitational constant 
and may self-interact through a potential A(<l>) ( the usual Brans-Dicke scalar field is 
(p = exp{—^) ). The function ci;($) plays the role of a coupling function between the 
dilaton and graviton. Each scalar-tensor theory is characterized by specific functional 
forms of a;(<l>) and A(<I>). A cosmological constant in the gravitational sector of the 
theory corresponds to the special case where A(<I') is a space-time constant. Among 
others we have the case of Brans-Dicke theory, where cc;($) is a space-time constant 
and A(<l') is absent [0]. It is known that inflationary solutions exist in a wide class 
of scalar-tensor cosmologies P, therefore these theories are relevant to the study of 
the very early Universe |^. Also, higher-order and higher-dimensional theories 
of gravity may be expressed in a scalar-tensor form after suitable field redefinitions. 
Brans-Dicke theory with cj = — 1 corresponds to a truncated version of the string 
effective action [^. 

Point symmetries associated with action (|l|) have been discussed previously within 

the context of the spatially isotropic Friedmann Universes P, ^ and Bianchi models 

T0[| . It was found that uj{^) and A($) must be related in a certain way if the field 



equations are to be symmetric. 

In this paper we consider the Wheeler-DeWitt equation (WDW) for the simplest 
anisotropic model, namely, the Bianchi type I cosmology. Some time ago the canonical 
formulation of the Brans-Dicke theory was considered by Toton and by Matzner et 
al. 11^]. For a recent review of WDW in different theories of gravitation see Ref. |jl2 



The line element for the class of spatially homogeneous space-times is given by 

ds"^ = -dt^ + habUJ^uj\ a,b = l, 2, 3, (2) 

where hab{t) is a function of cosmic time t and represents the metric on the surfaces 
of homogeneity and u"' are one-forms. These models have a topology R x G-s, where 



Gs represents a Lie group of isometries that acts transitively on the space-hke three- 
dimensional orbits |T^. The Lie algebra of G3 admits the structure constants Cbc = 
m'^'^edbc + (^"[feOc], where m"'' is a symmetric matrix, Qc = C^^ac and eabc = ^atc]- The 
Jacobi identity C'^bicC^de] = is only satisfied if m'^''afe = 0, so m"'^ must be transverse 
to a^ [l^. The model belongs to the Bianchi class A if a^ = and to the class B if 
Ofe 7^ 0. A basis may be found such that a^ = (a, 0, 0) and m"-'^ = diag [1^11,77122, ^^33], 
where ma take the values ±1 or 0. In the Bianchi class A, the Lie algebra is uniquely 
determined up to isomorphisms by the rank and signature of m°'^. The six possibilities 
are (0, 0, 0), (1, 0, 0), (1, —1, 0), (1, 1, 0), (1, 1, —1) and (1, 1, 1) and these correspond, 
respectively, to the Bianchi types I, II, VIq, VIIq, VIII and IX. Finally, the three- 
metric may be parametrized by hab{t) = e^"'-*^ fe^^^*-') , where e^" represents the 
effective spatial volume of the Universe and 

Pab = diag [p+ + V^p. , P+ - VSP^ , -2/5+] (3) 

is a traceless matrix that determines the anisotropy in the models. 

The configuration space Q for the Bianchi models derived from action (|I|) is there- 
fore four-dimensional and is spanned by {g„ = a,^,j3±}. The Lagrangiandensity 
L{qnAn) is defined hy A = J dtL{qn,qn), where a dot denotes differentiation with 
respect to cosmic time. It may be derived by substituting the metric @ into the 
action (P and integrating over the spatial variables. This procedure is unambiguous 



for the class A cosmologies and the action for these models simplifies to |jTO| 

A= f cite^"-* [6d<i> - 6a^ + 6/5^ + 6/?! + uj{<^)^^ - 2A($) + e-2"[/(/3±)] , (4) 

where 

U{(3±) = -e-^" (mabm'^' - -m^) (5) 

is the curvature potential, m = m°'a and indices are raised and lowered with h""^ 
and hab, respectively |T^. In the case of the type B models, a divergence may arise 
because the three-curvature contains a term proportional to aba^ ||16|. In view of this 
difficulty, we do not consider these models further. 

To proceed we take the following changes of variables. 



e^/^dr, x = a- $/2, y = J f-±^d<^, (6) 



and the action is now 



A 



dr 



6e3- {y'2 _ ^'2 ^ ^/^ ^ f^a_ _ ^(^)| ^ eXf/(^^) 



where 



A(y) := e*A($)/3. 
The above action means that the Lagrangian is 



L 



6e3- {y'2 _ 3,/2 ^ ^/2_ ^ ^.2_ _ ^^^^1 ^ gXf;^^^^ 



From here we can calculate the canonical momenta and the Hamiltonian 



(7) 
(8) 

(9) 



TT^ 



H 



Vle'^'x', -Ky = Ue-'^'y', n+ = 12e^^/5;, vr. = 12e''^/3l, 



-3x 



-^[-VT. 



n'y + vr^ + vrl + 144e^"A(y) - 24e^"?7(/3±)]. 



(10) 
fll) 



Then the Wheeler-DeWitt equation follows from the canonical quantization of if = 0, 
i.e., the canonical momenta in Eq.(p!0|) are converted into operators in the standard 
way, in general, owing to the ordering ambiguity we have 



I d f sd 

X 



therefore the WDW equation H'^{x,y) = for an arbitrary factor ordering, encoded 
in the B parameter, is 



52 _ j39^ _ Q2 -di-di- 24e^^t/(/3±) + 144e^^A(i/) vl>(a;, y, (3^ 



0. 



(13) 



In the following section we try to solve the Wheeler-DeWitt equation by separation of 
variables in the simplest of the homogeneous cosmologies, namely the Bianchi Type 
I case in which the potential U{(3±) vanishes. 



2 Exact solutions for Bianchi I 

We now consider the case of Bianchi type I cosmological model for which the potential 
U{j3±) vanishes identically. Furthermore we restrict ourselves to the case when A(y) = 
Ao/144 = constant, i.e. A(0) oc </>; this form of the potential has been used previously 



to obtain classical solutions in Bianclii type I vacuum cosmology [jT^ and for isotropic 



models with a barotropic fluid |T^, |T9|]; this potential is one of those that could produce 



infiation[^]. The choices that we have made simplify the WDW equation and allows 
us to obtain exact solutions by means of separation of variables, 

v^(x, y, /5±) = X{x)Y{y)F{(3+)G{(3.). (14) 

This implies that equation ([I3| ) gives the separeted equations 

X" - BX' + [kI + Aoe^"]X = 0, (15) 

Y" + k'^Y = 0, (16) 

F" + klF = 0, (17) 

G" + k^_G = 0, (18) 

where 

Kl = kl + k^_ + k\ (19) 

The solutions to these equations give the wavefunction: 



3 ^ ^ ^v 3 

c^e"'+^+ + Cfse~"'+'^+][cje'''-'^- + cge^'^-^-], (20) 



vl>(x,i/,/3±) = e^^/\rU^e'^) + c,Y,C^e'^)][c,e'y + c,e-^'y 



where 



-B - 4:t4 



6 ■ <"> 

and the q are constants. By superposition of these solutions, wavefunctions satisfying 
different boundary conditions can be obtained. In the following subsection we consider 
the case of wormholes. 



2.1 Wormhole solution 

Quantum wormholes can be regard as special class of solutions to the Wheeler-DeWitt 
equation with certain boundary conditions pO|- ^^^ wavefunction is exponentially 
damped for large spatial geometry, i.e., when a -^ oo, ii) the wavefunction is regular 
when the spatial geometry degenerates, i.e., the wavefunction does not oscillate when 
a -^ — oo. These conditions are known as the Hawking-Page regularity conditions 
(HP). In what follows we show, with a particular factor ordering, the explicit form of 
a wavefunction that satisfies the HP regularity conditions. 

The following wavefunction can be obtained by superposition of the solutions of 
the previous section or it can be substituted into the WDW equations to check that 
it is a particular solution 

where r is an arbitrary real parameter, the factor ordering B and the constants 
m, n,p, q satisfy the relations 

B = 3 + Ao/3, n^+p^ + q^ = 9, -{n^ + p^ + q^)m^ = Aq. (23) 

From these relations we see that for a wormhole we require a negative Aq. If from the 
last relation we take the negative root for m, it is easy to check that the wavefunction 
is exponentially damped for large spatial geometry, i.e., when a — i> oo (x — *> oo) and 
also that the wavefunction does not oscillate when a —>■ — oo( x —^ — oo). 

2.2 WKB solution 

Here we want to obtain the WKB wavefunction in the form 

^(^,e,/5±) = e'[^("'^'^+'^-)l (24) 

After substitution into the WDW equation the Hamilton- Jacobi equation results and 
separating variables it is straightforward to obtain the following solution. 



S = Pyy + V+I3+ + V-13- ± 



v/F + Aoe«- fcArcTanh(V^^^oZ:) 
kx 



(25) 



where k^ = P^ + p\ + P-- Once that we have solved the Hamilton- Jacob! equation 
it is possible to find the classical solution, we do not do that here because they were 
obtained by Banerjee et aLjlT^, solving the field equations. 



3 Another potential 



We consider now another potential function for which it is possible to obtain exact 
solutions to WDW equation. In Eq.(|^) we change the time r in the following way 
dr = e~^da, the action becomes, 



A = 6 da 



^"^{y"-x"+P'l + p"_}-e'^A{y) 



(26) 



here the prime means derivative with repect to a. We introduce a change of variables 
and an explicit potential, 



r^ = e'^Cosh(22/), ^ = e^^Sinh(2y), A{y) = AiCosh(2|/) + A2Sinh(2|/). (27) 

In the new time and variables, the action, Lagrangean, Hamiltonian can be calculated 
with the following results. 



A = 6 I da 

A 



2 (■2\/r/r 2 , nl 2x 



iC - v") + iv' - n{/3'+ + P'-) - Mv - A2^ 



-/2 ^/2\ I /^2 c'i\l n' 2 



:28i 



(29) 



L = 6[-(r' - v") + {v' - em + P'-) - AiV - A2e], 

H = ^[vrf - TT^ + (77^ - e){7i+' + nJ) + 24A,r] + 24A2e]. (30) 

The WDW equations in this case is 

'dl - dl + (r^2 - e){dl + dl) - 24Air/ - 2Ak,i] v&(r^,e,/?±) = 0. (31) 

Assuming separation of variables. 



vI/(r7,e,/5±) = E(r/)X(OP(/3+)M(/5_), 
the corresponding equations are 

E" + [k\rf + 24Air7 + k'^]E = 0, 



(32) 
(33) 



X" + [kii' - 24A2e + r]X = 0, (34) 



P" + klP = 0, (35) 



M" + r M = 0, (36) 

where k^,k- and k are arbitrary separation constants and k'j_ = k']^ + k'^_. The 
solutions are 



E = ci iFi(ai,l/2,2i) + C24/' iFi(ai + 1/2,3/2,^1), 

ai = ^M^ + (^f^/^+l],-l = ^/^(^ + ^)^ (37) 



X = C3 iFi(a2,l/2,22) + C4^2^' iFi(a2 + 1/2, 3/2,^2), 

«2 = ^[-A:^ + (^rV^+l],^2 = V^(e + ^)^ (38) 

P = c^e'^+P+ + c^e-'^+^+ , (39) 

M = C7e''^-^-+C8e-*'^-^^, (40) 

where iFi is the confluent hypergeometric function. More solutions can be obtained 
by superposition. 

3.1 WKB solution 

Again for the new potential we want to find the WKB wavefunction, 

^(^,e,/9±) = e*[^W+^(«)+(P+/3+)+(P-/3-)]. (41) 

After substitution into the WDW equation we have 

y^ - pW - 24Air7 - / = 0, (42) 



iy/-pie' + 24A2e-p^ 



where p\= p\ + p'^_. The solutions are 



(43) 



V = f^ + l] v'p2 + 24Air/ + p^2^2 + 



PX ^i^ "^ ^ % (44) 



-144Ai^+p2p^2\ jQg 



2p±^ 



VT 




p2 - 24 A2 r] + p^2 ^2 _^ 



(-144A22+p2p^2) log( '('^'^;^+""'") +2Vp^-24A2r7 + p^^r7^) 



2p±=^ 



.(45) 



4 Final remarks 



In this work we have obtained exact solutions for the WDW equation using the general 
scalar tensor theory of gravitation with arbitary coupling function uj{<^) and two 
specific selfinteracting potentials A($). The WKB wavefunctions were obatained for 
both cases. The solutions here obtained are for arbitrary large anisotropies. We hope 
that these solutions should prove to be useful in the study of the issues of quantum 
cosmology. In the past quantum cosmology using general relativity in homogeneous 
spacetimes with small anisotropies and a scalar field were considered by Lukash and 
Schmidt |2l|], and Amsterdamski [^. Lidsey pSf has considered the wavefunctions 
in a highly anisotropic cosmologies with a massless minimally coupled scalar field. 
More recently Bachmann and Schmidt [^] have considered the case of arbitrary 



anisotropies in Bianchi I quantum cosmology. 
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